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Abstract. We introduce stochastic models for the transport of heat in systems described 
by local coUisional dynamics. The dynamics consists of tracer particles moving through 
an array of hot scatterers describing the effect of heat baths at fixed temperatures. Those 
models have the structure of Markov renewal processes. We study their ergodic properties 
in details and provide a useful formula for the cumulant generating function of the time 
integrated energy current. We observe that out of thermal equilibrium, the generating 
function is not analytic. When the set of temperatures of the scatterers is fixed by the 
condition that in average no energy is exchanged between the scatterers and the system, 
different behaviours may arise. When the tracer particles are allowed to travel freely 
through the whole array of scatterers, the temperature profile is linear. If the particles are 
locked in between scatterers, the temperature profile becomes nonlinear. In both cases, 
the thermal conductivity is interpreted as a frequency of collision between tracers and 
scatterers. 



1.1. Lattice Hamiltonian dynamics. 

Smooth interactions. The study of conduction of thermal energy by Hamiltonian lattice 
dynamics has recently known a great deal of activity and numerical as well as analytical 
results have accumulated. For one-dimensional systems, general lattice Hamiltonian dy- 
namics may be defined in the following way. Consider N particles of unit mass located on 
a one-dimensional lattice with local positions and momenta (q, p) = {(qj, pj)}^^.^^, with 

qj,Pi G M'^. The Hamiltonian H takes the form, 

^ r„2 1 



where V represents the interaction with an external substrate and U a nearest-neighbor 
interaction. Typically, the pinning potentials V and U are smooth and confining, i.e. 
they grow to infinity when the norm of their argument goes to infinity. Much effort has 
been devoted to the study of those systems in the case where the dynamics is a small 
perturbation of a completely integrable one, namely the case of chains of weakly anharmonic 
oscillators. For particles moving in a one-dimensional space, those are described for instance 
by potentials of the form. 



with A small . The situation may be summarized as follows fL\ \6[ [T^ I19j. the dynamics is 
properly described by a Peierls-Boltzmann equation for phonons and under the assumption 
that this equation holds, the conductivity is finite and may be computed as a function 
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of temperature and microscopic interactions. The temperature dependence is of the form 

i/r2. 

CoUisional dynamics. Another type of lattice Hamiltonian dynamics arises in systems where 
the local dynamics is given by a billiard dynamics ^71ll H [T2 l ll6j. Physically, dynamics of this 
type can model aerogels, namely gels from which one has removed the liquid components 
and replaced them by molecules of gas. For instance, one may take in (jl.ip particles moving 
in one dimension, i.e d = 1, and take a sequence of interactions and Uk, 

2k 

Vk{x) = h {D , Uk{x) = h (^) , fk{x) 



2k 



In the limit k ^ oo, one obtains, 



\ Oif |x| <6 ~ \ Oif |x| <a ^ ' 

and the dynamics is described by a sequence of "collisions" between nearest neighbors. 
When the difference of positions of two neighboring particles reach the parameter a, they 
exchange their velocities and thus their kinetic energies. This model, originally introduced 
in [16], was dubbed in [12] the complete exchange model. One may also simply consider 
particles moving in square cells located on a one-dimensional lattice. While remaining 
confined at all times, the particles collide with their nearest-neighbors through holes in 
the cells walls. The first type of such models, in which the local dynamics is described by 
semi-dispersing billiards was introduced in [7j and its thermal transport properties were 
studied in [TT]. General collisional models of this type have been introduced in [12] and 
may be described as follows. Formally, the dynamics is described by a Hamiltonian of 
the form (II. Ih where interaction potentials is equal to zero inside a region fif/ C with 
smooth boundary A of dimension d—1, and equal to infinity outside. Likewise, the pinning 
potential V is assumed to be zero inside a bounded region and infinity outside, implying 
that the motion of a single particle remains confined for all times. The regions 0[/ and 
being specified, the dynamics is equivalent to a billiard in high dimension. 

Physical observables. Dynamics described by smooth interactions and collisional dynamics 
have several qualitative differences and similarities that we briefiy point out. In both cases 
the interactions between components occur only between nearest neighbors on a lattice and 
thus the evolution of the local energy may be written as, 

En{t) - En{0) = Jn-l^n{[0,t]) - J„_„+i ( [0, t] ) . (1.3) 

In the case of smooth interaction potentials, the time integrated energy current between 
sites n and n + 1 takes the form 

J„^„+i([0,t]) = / -(p„(s) + p„+i(s)) • VC/(q„(s) - q„+i(s))ds 
Jo ^ 

whereas in the case of collisional dynamics, 

J„_„+i([0,t]) = \Y. {vn^Slf-vU^{Slf\, (1.4) 

0<k<Nt 

where the component of the vector p„ in the direction of the unit vector n = H^i — 
(7i+i||~^((7j — (7j+i) at the time of collision is denoted by Pn = Pn ■ n. Nt counts the 
number of collision up to time t and {S^)k is the sequence of collision times. Note that 
in the complete exchange model, p:}^ = pn and the time integrated current between two 
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neighbors is simply the sum of ah kinetic energy exchanges between the particles. Assume 
now that such systems are thermalized at different temperatures at their boundaries. In 
order to understand the transfer of energy from one side to the other, one is interested 
in the ergodic behavior of the current and in computing lim^^+oo i~^«^ra^n+i([0, t]) which 
gives the average current of energy in the stationary state. Because of the special form of 
the time-integrated current <\lAh . a natural guess to make is that local equilibrium settles 
in and to assume that this limit is given by 



where = linit^oo t Nf is the frequency of collisions between neighbors under local equilib- 
rium conditions. r„ = ^{Pn) is the average kinetic energy of the particles. The conductivity 
is thus identified to the frequency of collisions. Typically, the collisions occur when the par- 
ticles get near the boundaries of their cell and thus the frequency of collisions is roughly 
proportional to the average number of visits to the boundaries per unit time. Because the 
particle travels freely within its cell, this is proportional to \/T^. Once this relation is taken 
for granted, then the temperature profile may be computed because, by conservation of 
energy, the current must be constant throughout the system, 



This is equivalent to a finite difference equation for the set of temperatures of the local 
equilibrium distribution. Numerical studies show that the identification of the conductivity 
with the frequency of collisions holds true to a very high degree of accuracy in a wide class 
of collisional dynamics, when the individual particles collide rarely. In particular, it was 
shown [12J that it does not depend on the detailed chaotic properties of the local dynamics 
as it was originally assumed [11]. 

Although such a simple way to compute heuristically the heat conductivity does not exist 
in weakly anharmonic systems, there are similarities between the two dynamics considered 
in a weakly interacting regime and written in proper coordinates. In both linearized 
Boltzmann equation describes properly the thermal properties of the systems. Depending 
on the context, each mode or particle behaves as if it was coupled to an ideal stochastic heat 
bath and the intensity of the coupling depends on the microscopic interactions and yields 
the thermal conductivity. This quantity is also identified with the frequency of collisions 
between the components (phonons or particles) in equilibrium. More precisely in the case 
of collisional dynamics, each particle moves freely within its own cell and interact with 
its neighbors as if those were part of infinite thermal bath with fixed temperature. The 
Boltzmann type approach is successful in computing theoretically the conductivity from 
microscopic interactions and local temperature as was checked in numerical simulations 



1.2. Models and results. An important feature of the collisional models is that the evo- 
lution of energy occurs at discrete (collision) times and amounts to an exchange of kinetic 
energy between neighbors. The length of the interval of time between two successive colli- 
sions depends itself on the kinetic energy of the particle, which fixes the VT dependance of 
the collision frequency. Our idea is to build and analyze models which are stochastic from 
the start and share the general structure described by the Boltzmann description sketched 
above. 




(1.5) 



lim -J(„_i)^„([0,t])= lim - J„_^„,+i([0, t]). 

->+CXD t t^ + OO t 



mm- 
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Thus, we want to consider dynamics which consists of a mixture of integrable Hamilton- 
ian dynamics and colhsions with stochastic heat baths. The models are made of scatterers 
described as heat baths and tracer particles transferring energy between those "hot" scat- 
terers. The tracers move in a one-dimensional interval in which the scatterers are located 
on a lattice. The motion of a tracer is ballistic except when it encounters a scatterer. 
At that point its velocity is randomly updated according to a law which depends on the 
temperature of the scatterer. The temperature of the scatterers is fixed by the condition 
that in the stationary state, no energy is exchanged between the scatterers and the tracers. 
In spirit, this is similar to the so-called self-consistent chain of (an)harmonic oscillators 
[U m [5], but in our case it may be naturally interpreted as a condition ensuring that the 
energy transfer per unit time between scatterers is constant throughout the system. 

Relation to previous works. Geometrically, the systems we study are analogous to the ones 
introduced and studied in [9l \T0\ I15j . The scatterers of our models are similar to the 
energy storing devices of those models. However, in our case, the dynamics is stochastic 
from the start and the action of the scatterer models the one of a very large system. 
As in those models we distinguish two types of dynamics. Depending on whether the 
particles are confined or not between the scatterers, we derive a non-linear or linear profile 
of temperatures for the scatterers. Those two types of behaviours described by wandering 
or confined tracers seem to be universal in systems described by a collisional dynamics at 
a microscopic level [U [101 El dH [I2l [16] . Among the collisional dynamics described above, 
two typical examples are given in [16] and [111 [T2] . In [16] , a detailed numerical analysis of 
the complete exchange model is provided. Fourier's law holds and the temperature profile 
is linear. In our framework this may be understood as an instance of wandering tracer 
dynamics. Indeed, in the complete exchange dynamics, it is not only the energy that is 
exchanged between neighbors but also the momenta of the particles. Thus, the dynamics is 
more similar to the one of particles traveling through the whole system. In |1H I12j. models 
where the particles are confined and exchange only a fraction of their energy are considered. 
They display a temperature profile identical to the ones of the confined tracers. 

Organisation of the paper. As stochastic processes, our models are naturally described by 
Markov renewal processes and in section [2] we consider the simplest dynamics of the type we 
want to study. We define the stochastic process generated by the free motion of a particle in 
a box. The particle collides with the walls of the box and its velocity is randomly updated 
according to some fixed probability law. We compute the unique invariant measure for 
this process. This justifies rigorously the updating rule used in computer simulations |20j . 
In section [3l we define the tracers and scatterers models in full generality and provide an 
explicit formula for the stationary measure out of equilibrium. In section U we describe 
and analyze the properties of the wandering tracer model in its simplest version, i.e. when 
a tracer encounters a scatterer it is deterministically transmitted on the other side of the 
scatter. We find that when the temperatures are fixed and the temperatures in the bulk 
chosen such that the transfer of energy is constant throughout the system, then the tem- 
perature profile of the scatterers is linear. The identity between thermal conductivity and 
frequency of collisions between tracers and scatterers appears as a natural consequence of 
the renewal theorem for Markov renewal processes. Next, we study the cumulant generating 
function of the time-integrated current of energy. We give a rather explicit formula allowing 
to compute derivatives of any order. A striking feature is the lack of analyticity of the the 
generating function. The origin of this phenomenon may traced back to the presence of 
particles with arbitrarily low speed. Nevertheless, we are able to show the validity of the 
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Green-Kubo formula for the conductivity. The final section [S] is devoted to the analysis of 
the dynamics of confined tracers. In that case, we find that the condition that there are no 
exchange of energy between the scatterers and the particles imposes a non-linear profile of 
temperature. 



2. Basic process. 

We first describe in detail a model which will be the elementary building block of the 
models we intend to study. We consider a single particle moving in the interval [0, 1[ with 
a positive velocity. When the particle reaches 1, it is absorbed and re-emitted in with a 
random positive velocity. 

To be more precise, we consider an i.i.d. sequence (ui)i=i,2,... such that Vi > a.s. for all 
i. The particle is re-emitted from with speed Vi after its i-th collision with the wall at 1. 
The time to reach 1 again is then Tj := 1/vi. 

Suppose the particle starts at time t = at position qq G [0, 1[ with velocity po > 0. The 
time of the first collision with the wall at 1 is 

So = So{qo,Po) ■= - — —, 
Po 

and the time of the n-th collision is 

Sn ■■= So + Ti -\ h r„, n > 1. 

We define now a stochastic process {qt,Pt)t>o with values in [0, l[xM_|_ 

{qo + pot,p) if t < So, 

{qt,pt) = F{q,p,t,{Tn)n>l) ■ = 

if Sn-l<t<Sn, n>l. 

We denote by Bh the set of all bounded Borel / : [0, l[xR+ ^ E and we set 

Ptf{qo,Po) :=nf{qt,Pt))=miFiQ.P,t,{Tn)n>i)), {qo,Po) G [0,l[xR+. 
Then it is not difficult to prove that 

Proposition 2.1. The process {qt,Pt)t>o Markov and {Pt)t>o has the semigroup property: 
Pt+s = PtPs,t,s>0. 

2.1. The invariant measure. We assume now that 

fi := E{Ti) = E{l/vi) < +00, 

and that the distribution of Tj is non-lattice, i.e. there is no (5 > such that P(Tj G (5N) = 1. 
This assumption is not necessary but it simplifies the presentation; for the applications we 
have in mind, the distribution of Tj has a density and is therefore always non-lattice. 

We denote the law of Tj by ip{dT) and the law ofvi = I/tj by <l){du). The law of ti + - • ■+Tn 
is denoted as usual by the n-fold convolution ip^* . 

Proposition 2.2. The only invariant measure on [0, l[xM-|- of the process {qt,Pt,t ^ 0) is 
given by ^{dq,dp) = dq (j){dp)/ fip. 
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Proof. For any bounded Borel function / : [0, l[xM_|_ i— > R 
Ptf{qo,Po) = 

= l(t<So)/(^0+P0i,P0) + l(t>So) \'^{Sr.-l<t<Sn) f ^]] 

n=l ^ \ -Tn TnJ J 

= l(t<s„) /('Zo + Pot,Po) + l(t>s„) E / ^("-'^(ds) / ^{dr) f , - 



l(t<s„) f{qo + Pot,Po) + l(t>s„) U{ds) J ^{dr) f i 



T T 



where we recall that Sn = So{q,p) + ti + ■ ■ ■ + Tn and we set ip (ds) = 6o{ds) and 

/b °° fb °° pb 

U{ds) = i^*^''~^\ds) = 5o{ds) + < a < 6. 

n=l n=l 

The renewal measure U{ds) gives the average number of collisions in the time interval ds. 
We obtain 

Ptf{qo,Po) = i{t<So)fiQo+Pot,Po) + lit>So) ^{dr) Utids)fl-,-j (2.1) 

where Ut{[a,b]) = U{[t — b,t — a]) for < a < b <t. By Blackwell's renewal theorem (see 
e.g. [U Theorem V.4.3]), U{[t — b,t — a]) converges to {b — a)/ fi as t ^ +oo. Therefore, the 
last expression converges to 



= f f{x,u)— — ^- dx =: [ f{x,u)^{dx,du). 

J[0,lfxM+ J[0,lfxR+ 



/i \ T 



'[0,l[xR+ J[0,l[xR4 
In other words we have for all bounded Borel / : [0, l[xl 



Jim Ptf{qo,Po)= [ fd^, V ((7o,Po) G [0,l[xM+. 

Since {Pt)t>o is a semigroup, we obtain 

/ fd'y= lim Psf{qo,Po)= lim Pt+sfiqo,Po) = lim PsPtf{qQ,Po)= / ^t/c?7, 

J s— >+oo S— ^+00 s— >+oo J 

i.e. 7 is invariant for {Pt)t>o- This convergence result for all initial conditions {qo,Po) £ 
[0, l[xR+ implies that j{dq,dp) is the only invariant measure of the process. □ 

We remark that the invariant measure 7 satisfies 



fdj = E^f{U,T-') -j, (2.2) 
where {U,t) is independent, U is uniform on [0, 1] and r has same law as ri. 
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2.2. The main example. The process we study here will be used to build dynamics where 
the updating of the velocities simulate the action of a thermal bath acting through a wall. 
Assume that in our set-up one wishes to obtain a process whose invariant measure jidq, dp) 
is given by 

l{dq,dp) = l[o,i](g) l]o,+oo[(p) V — e"^^ dqdp, 

V TT 

where f3 > 0. Then from the above results we see that we must choose the updating 
distribution 



Mp) = hp>o)P\Me-P'^ ^ = l(p>o)/?pe-'^^ . (2.3) 

The distribution of the interarrival time becomes accordingly 

ipp{dT) := l(^>o) ^ exp dr, (2.4) 

where we recall that (fip and ilip are related by 

f{v)ct)p{v)dv = j ^p{dT). 

2.3. A renewal theory viewpoint. We recall now some classical results of renewal theory 
and show their relevance in our setting. 

We consider again the i.i.d. sequence (t„)„>i that we used for the construction of the 
process {qt,Pt)- We recall that r„ > a.s., E(r„) := /u < +00 and the distribution of Tj is 
non-lattice. A (delayed) renewal process associated with (r„)„>i is a sequence 

"S'o > 0, Sn ■■= Sn-l + Tn, n > 1, 

where S'o is independent of {Tn)n>i- We define 

00 

m := ^ l(s„<t) = #{"- = 0,1, . . . : Sn <t} = inf{n : Sn > t}. 

n=0 

We assume for simplicity that 5o > and Sq is non-random. Then we define the processes 

Bt := Sn, -t, t>0, 

Ao + t, t<So, 

where ^0 ^ is non-random and Bq = Sq > hy assumption. Then it is easy to see that 
in our setting, for {qo,Po) G [0, l[x]0, -|-cx)[ 

^0 = —, Bo = So = So{qo,po) ^ qt = -r-^^, Pt = . "| „ , V t > 0. 
Po At + Bt At + Bt 

This follows, in particular, from 



„2 



1 



Po, y t < So, 



At + Bt 



Tnt^ y t> So, 

Then we recall the following result (see e.g. O Chapter V]) 
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Proposition 2.3. The processes {At)t>o, (-Bt)i>o ^.i^'d {At,Bt)t>o are strong Markov. Any 
invariant probability measure of {At)t>o is also invariant for {Bt)t>o and must coincide 
with xP(ti G dx)/iJ,. The only invariant probability measure v for the process {A^, Bt)t^Q 
is given by 

J fdv=E(^f{UT,{i-u)T) 

where {U,t) is independent, U is uniform on [0, 1] and r has same law as ti. 

Since {qt,Pt) is given by an injective function computed at {At, Bt), then {qt,Pt) is also strong 
Markov; moreover, an invariant probability measure 'y{dq, dp) for {qt,Pt) is necessarily given 

which is the same formula found in ()2.2p above. 

2.4. Markov renewal processes. A natural generalization of a renewal process is pro- 
vided by a Markov renewal process (see |2j), which is a pair of stochastic processes {Xn, Tn+i)n 
such that 

(1) {Xn)n>o is a Markov chain in some finite state space E with transition density 

{Qij)i,jeE 

(2) conditionally on 7i = ct((X„)„>o), (t„+i)„>o is an independent sequence of positive 
random variables such that for all n > 1 

P(r„ <t\n)=F{Tn<t\Xn, Xn+l) = Gij{t) 

on the event {Xn = i,Xn+i = j}, where {Gij)ij(zE is a family of distribution 
functions on ]0, +cxd[. 

The sequence of interarrival times (rn+i)n>o is in general not i.i.d. However, given the 
Markov chain (X„)„>o, {Tn+i)n>o is independent and the law of Tn+i depends on X^ and 
This kind of process will be the main mathematical tool in the concrete models we 
analyze in the following sections. 

3. Tracers and hot scatterers 

We consider a gas of M non-interacting tracer particles moving through a one-dimensional 
lattice of scatterers Wn, n = 1, . . . ,N . In between the scatterers, the tracers move with 
constant speed in the boxes /«, 

I„ = [n-l,n], n = l,...,iV (3.1) 

When a tracer encounters the scatterer n, it is absorbed and re-emitted on either side ac- 
cording to a certain probability distribution with a random velocity p distributed according 
to the law 

4>tSp)=P^Pne-f'-"^, (3.2) 

see ()2.3p . The state space describing the motion of the particle is thus the cartesian product 
of the positions space / = [0,A^] and velocity space R* = M\{0}, = / x M*. At the 
extremities of the system the sign of the velocity is reversed but the particle bounces back 
with a random velocity distributed according to the same law with parameter (3l and (Sr. 
This dynamics defines a Markov renewal process which we describe now more formally 
using the same notations as above. 
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For notational simplicity, we define the dynamics for a single tracer particle, the extension 
to M particles is straightforward. The particle moves in the interval I := [0,N], which is 
split into N subintervals of equal length: In := [n — l,n], n = 1, . . . , N . At time t = 0, 
the particle starts at position qo £]0,N[ with speed po G M* := K\{0} and we define 
(no, o"o) G E as follows: no := [q + (sign(po) + l)/2j G {0, 1, . . . , N}, where [-J denotes the 
integer part; in other words no is such that 

no - 1 < go < -^0 if sign(po) = +1, 
no<qo<no + l if sign(po) = -1- 

and denotes the first scatterer the particle hits. We define moreover 

f sign(po), if no G {!,... ,iV - 1}, 
\ -sign(po), if noG{0,iV}. 

In other words, if the first scatterer the particle hits is at the boundary {0, N} of the system, 
then the particle will be reflected at the first hitting. 

We suppose that the sequence = {nk,crk) of scatterers visited by the particle and 
signs of the velocity is a Markov chain on 

E:={{n,a), n=l,...,N-l, a = ±1} U {(0, +1), (iV, -1)}, 

with initial state (no,o"o) and with an irreducible probability transition matrix on E such 
that 

( 1 if (n,f7) = (1,-1) and (n',cr') = (0,+l) 
(l{n,aUn',a') = < 1 if i^i, a) = (N - 1, +1) and (n', a') = (N, -1) 
[O if nG {l,...,iV-l} and n' - n / a'. 

The first two conditions mean that the tracer is always reflected at n = and n = N. The 
last condition means that the new sign a' gives the next scatterer n' visited: if a' = +1 then 
n' = n + 1, if fj' = — 1 then n' = n — 1. The irreducibility assumption gives the existence of 
a unique invariant probability measure that we call {i^a)a£E- 

We now define the time the particle takes between two subsequent visits to the scatterers. 
Conditionally on Ti. = o-{{Xk)k>o), the sequence {Tk)k>i is independent with distribution 
defined by 

Pn ( I3n 



P(t,, GdT\n)= ¥{Tk G dr I Xk-i) = ^ exp l(,>o) dr =: Vn(fir) (3.3) 



on the event {Xk^i = (n, o")}, where /3o, • • • , Pn G M+, see ()2.4p . 

We consider now the Markov chain {Xk)k>o with initial state Xq = {n,a) and the 
associated sequence {Tk)k>i- The time of the first collision with a wall is 

'^'0 = ^o{qo,Po) ■■= > 0, 

Po 

and the time of the A;-th collision with one of the scatterers is 

Sk ■■= So + Ti-\ hTfc, A; > 1. 

Before time 5*0, the particle moves with uniform velocity p. Between time Sk-i and time 
Sk, the particle moves with uniform velocity — and (Sk)k>o is the sequence of times when 



G {0, . . . , A^}. In particular we define the sequence of incoming velocity Vk at time S, 



k 



vo ■=Po, Vk :=—, k> 1. (3.4) 

Tk 
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We define precisely the stochastic process {qt,Pt)t>o with values in [0, A^] x M* 

{qo+Pot,Po) if t < So, 

{qt,Pt) ■= { . s (3-5) 

^ (nk-i + f (t - Sk-i), f ) if Sk-i <t<Sk, k>l, 

and we use the notation Xj^. = {rik^ak), ^A;-i = (^^fe-ii Then, in analogy with 

Propositions I2.l l and 12.2] we have the result 

Proposition 3.1. The process {qt,Pt)t>o Markov and its only invariant measure on 
[0, 1] X M* is given by 



lidq, dp) = —^ tiM) ^ ( ^{n-l,a) 9(n-l,^),{n,+l) l(p>0) Pn-l e"^""! V + 

(3.6) 

+ T^{n,a) q{n,a),(n-l-l) l(p<0) Pn C"^"^ J dq dp 



'-^ n=l a=±l 



where Zn = x/f E(n,a)Gi? VK- 



Proof. The proof is very similar to that of Proposition 12.21 The renewal measures we have 
to consider is 



Ua,a'it) = J2^a {Sk < t, = a') , 



k=l 

where under Pq, the process X^ has the law described above with Xq = a = (i, a) a.s. Then 
the analog of formula (j2.1|) is 

Ptf{qo,Po) = l(t<So) fiQo+Pot,Po) 

/ qa',a" i^n'idr) / Ui^^,{ds) / n' H , — 

where we use the notation a = (n, a), a' = {n' , a'), a" = {n" , a") and the measure ^fids) 
are defined by 

UiA[a,b]) = Uc,At-b)-UaAt-a), 0<a<b<t, a,(5^E. 

Since the analog of Blackwell's Theorem holds also for Markov renewal processes, see e.g. 
[21 Theorem VII.4.3], we obtain 

hm [/* ([a,6]) = (6-a)^ 

t— >+oo ^ 

where v is the unique probability invariant measure on E of the Markov chain {Xk)k>Q and 
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Therefore 

lim Ptfiqo,Po)= / Qa' ,a" i^n' {dr) ds^fin'^ , — 

^ /•I r+co / cr" \ 

= = '^a'qa',a" dx T ipn' (dr) f ( n' + a" X , — ) 

^a',a"GE Jo Jo V 

f{q,p)l{dq,dp). 

'[0,l]xK 

Arguing as in the proof of Proposition 12. H we conclude. □ 

Remark 3.2. Notice that the invariant measure is always explicit, although in general the 
process {qt,Pt)t>o is not reversible. 

4. Wandering tracers 

4.1. Generalities and physical observables. In the first model that we study, when a 
tracer reaches a scatterer n G {1, . . . , — 1}, it is absorbed on one side and re-emitted 
on the other side with a random velocity distributed according to a law determined by the 
temperature of the scatterers. The sign of the velocity changes when and only when the 
tracer reaches the scatterers or N. The transition matrix of the underlying Markov chain 
is 

' 1 if n' = n + a ^ {0,N} and a = a' 
1 if (n, a) = (1, -1) and (n', a') = (0, +1) 
1 if (n,f7) = (A- 1,+1) and (n',cj') = (iV,-l) 
otherwise. 



q{n,a),{n',a') 



The associated invariant probability measure given is the uniform distribution on E. In 
fact, in this case the Markov chain moves deterministically as follows: A^ = {n/., Gk) where 

nfc = /(|no + o-oA;|mod2A), f{i):=N-\N-i\, i = 0, ...,2A 

(4.1) 

Cjfc = do (-1) L(n+-o(no-JV)+7V)/7VJ ^ ^ > i_ 

In particular, we have the following periodicity 

(Afc)fe>o = (Afc+2iv)/c>o under IP(no,(jo) 

(4.2) 

where = denotes equality in distribution. Proposition 13. II becomes 

Proposition 4.1. The process {qt,Pt)t>o is Markov and its only invariant measure on 
[0, 1] X R* is given by 

7idq,dp) = ^Y1 '^^r.iq) (l(p>0) Pn-i e-^-i'^ + l(p<o) Pn e'^-'^') dqdp (4.3) 

n=l ^ ' 

where Zn = Y.n=i ( V f^n-i + 



The first useful result is the computation of the asymptotic frequency of collision of a tracer 
with a fixed scatterer. 
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Proposition 4.2. For n € {0, . . . , A'"} set 4)^,0 '■= inf{£ > : = n}, 



El(-V.^*)' ^"-E2 1(5. *>0. (4.4) 



t ■— 

k=l k=l 



Then for any initial condition {qo,Po), ^(qo-po)''^-^- 

lim lim -^ = ^, if n G {1, . . . , iV - 1}, (4.5) 

iVf 1 

lim ^ = — , if nG {0,A^}. (4.6) 

t^+co t Z]\f 

Proof. Let n £ {1, . . . ,N - 1}. By ([12]), the sequence (5'</,„_2(fc+i) ~ 'S'0„.2Jfc>o is i-i-d. and 
therefore we have by the renewal theorem 



hm 



Since \Np — NJ^\ < 1, we conclude. □ 

We next identify the physical quantities of interest. The energy exchanged between the 
scatterer n and a particle during a time interval [0, t] is given by 



2 / ^ \~K~I-L ~K,) -'■(nfc=?i)) 

fc>0:Sfc<t 

recall that, by ()3.4p and (j3.5p . t;^ and f/t+i are respectively the incoming and the outcoming 
velocity at time 5^. The total entropy flow due to the exchange of energy between the 
scatterers and a particle is given by 

5n([0,i]) := -MM), 5([0,t]) := 5^5„([0,t]). (4.7) 

n=0 

The energy exchanged between scatterers n and (n + 1) during a time interval [0, t] is 

J„^„+l([0,t]) := - ^ <7fe_i=l) - l(nfe_i=n+l, afe_i=-l)) ■ 

k>\:Sk<t 

We define the energy flow per unit time in the stationary state by 

En:= lim y^„([0,t]). (4.8) 

Similarly, the entropy flow per unit time is given by 

1 ^ 
5„:= lim -5„([0,t]), 5:=V5„ (4.9) 

t^+oo t ^ ' 

n=0 

and the current of energy between scatterers Wn and tf^+i is given by the transfer of energy 
per unit time, 

Jn--= lim 7 J„^n+i([0,t]). (4.10) 
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Proposition 4.3. The limits in (|4.8p . (j4.9p and (j4.1Up exist ^(q^pg) a-S- and for all n = 

1,...,N-1, 

£n = ^ and £o = — , £n = , (4.11) 

Jn = ^^^, 5^ lg(7kz^>0. (4.12) 

n=0 

Proof. Setting Yk := {Xi+2Nk,i = 0, . . . , 2iV - 1) G E^^ then by (jM]) and the Markov 
property we have that (Yk)k>o forms an i.i.d. sequence. Using the notation (j4.4p we can 
write 



2 _ 2 

fc>0:5fe<i fc=0 



and define 

fc 2 



2 _ 2 



Eni[0,t])= X; Kfc+i + e^O- 

fc=0 



Notice that is the exchange of energy between the scatterer n and the particle at the (/c + 
l)th passage of the particle by the scatterer n. Since {Yk)k>o is i.i.d., then {e2f._^_l + e^k)k>o 
is also i.i.d. and we obtain by (j4.5p and by the law of large numbers for n G {1, . . . , — 1} 

1 - iV" 2 - 1 

hm -i?„([0,t])= hm ^ ii;„([0,t]) = — E(,„,p„)(e^ + eS). 

Now, over a period of 2iV transitions of {Xk)k>o, each scatterer n G {1, . . . , A^ — 1} is visited 
twice, once coming from the right and once from the left. Therefore 

E{go,Po)(e? + eo) = 2r„ - Tn-i - r„+i, V n G {1, . . . , - 1}. 

Notice now that |£^„([0,t]) - En{[0,t])\ < |e^„| =: Wt- It can be seen that Wt/t and 
therefore we obtain the first relation of (|4.1ip . 
If now n G {0,iV}, by (gSI) 

1 iV," 1 11 

hm -En{[0,t])= lira ^ ii;„([0, t]) = - — E(,„,p„) (e^) . 



In particular 



2 ^(-/o.Po) (^o) = ?o - Ti, - (e^) = Tjv - Ta^.i, 



and the proof of (|4.1ip is complete. The proof of the first relation in (|4.12p is similar. The 
second relation follows from ()4.1ip and a summation by parts: 



N-l f 1 1 \ 

n=0 n=0 



\2 



□ 
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4.2. Self-consistency, Fourier's law and temperature profile. From (j4.1ip and (|4.12p . 

we have the obvious result 

Proposition 4.4. (Self-consistency condition) The only collection {Tn)n=o,...,N such 
that 

£n = Q, n = l,...,N -I 

with To = Tl and T^r = Tr is 

Tt 

Tn = TL + -{TR-TL), n = 0,...,N. (4.13) 
The entropy flow per unit time S is equal to if and only if Tl = Tr . 

Note that the condition on the exchange of energy is imposed only for the scatterers. In 
contrast, when Tl ^ Tr the tracer will always exchange energy with the boundary walls. 
Let us consider now Mjv (to be fixed) non-interacting tracer particles described by their 
momenta and positions {p,q) = {pi,qi)i<i<MN moving through the array of scatterers. 
As the motions of the tracers are independent, the generalization is straightforward. The 
corresponding stationary measure is given by 



7^^(dg, dp) = n lidqu dp^) (4.14) 



i=l 

and the total average current between scatterers Wn and Wn+i is the sum of the contribution 
of each particle in ()4.12p 

jM. = Mn^^!^^^. (4.15) 

The total rate of energy exchanged between the scatterer Wn and the tracers in the station- 
ary state is given by 

g^^N = Mn ~ \' ~ . (4.16) 

Thus, if the the self-consistency condition is imposed and the temperatures of the scatterers 
is given by ()4.13p . then, by (j4.15p . one has 

'^n''' = -^iTR-TL). (4.17) 

The local conductivity is defined as the ratio of the average current of energy to the local 
temperature gradient, namely, 

K„= lim ^"^"^ = hm — . (4.18) 

N^oo Tn — Tn+l N-*oo Zn 

If the temperature profile is given by (I4.13p . then one may compute the explicit asymptotic 
behavior of Zn in the large limit 



lim ^ = ^2^ r , = (4.19) 

N^oc N Jo {Tl + x{Tr-Tl)Y^ Tl+T'l 

Thus for a number of tracers Mn = o{N), we have Hn = 0. This is because when the 
number of scatterers increases, the proportion of time that a given tracer spends carrying 
energy from scatterer Wn to Wn+i goes to zero simply because the tracer must go back and 
forth in a larger and larger system. However, we see that if we take as many tracer particles 
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as scatterers, namely M/v = N, then Fourier's law holds, i.e the conductivity is finite. Its 
value is given by 

1 i_ 

2 _|_ 2 

- (4.20) 



2V27r 

Notice that does not depend on n. In particular, \i = Tji = T 

4.3. Cumulant generating function and the Gallavotti-Cohen symmetry relation. 

We denote B = (/3o, • • • , /3iv) G M;'^^^. We do not need in this section to assume that 
Tj := j3~^ satisfy (|4.13p . Fix n G {0, . . . , — 1}. We are going to compute and study the 
properties of the cumulant generating function, 

fn{\B):= lim ^ logE(exp(-AJ„([0,i]))), VA G ] - (4.22) 

We define A„(q!) := (l(j=„, ^=i) - l(j=„+i, <7=-i)) for a = {i,a) G E. For A g] - f3n,Pn+i[ 
and e > 0, we define 

r+oo 2 
Cn{a,X,e) := Pi v 6-^-'-'^'+^^"''"^^ dv, a = {i,a)eE. (4.23) 
Jo 

and the function Fn, which is crucial in the computation of fn{X,B), 

Fn{\e,B) := II C„(a,A,e). 

Fn will be identified with the spectral radius of some matrix. 

We anticipate a striking feature of our model: for /?„ 7^ (3n+i, the cumulant generating 
function fn{-,B) is not analytic around A = 0. In particular, fn{-,B) > for A in a left 
neighborhood of 0, while fn{-,B) = in a right neighborhood of 0. See Remark 14.61 and 
section 133] below for further discussions. 



Proposition 4.5. < then VA g] - /?.„, 0[U]/3„+i - /3n,Pn+i[, fn{\B) is given 

by the unique solution €0 > to the equation 

Fn{X,€o,B) = 1. 

If X £ [0,Pn+i — Pn], then fn{X,B) = 0. The function fn{--,B) is convex and continuous 
over] — (3n, Pn+i[ o-nd satisfies the Gallavotti-Cohen symmetry relation 

fn{X,B) = fnWn+l " /3„ " A,^). (4.24) 

Proof. We call a family of cr-finite measures F^^a' (d-T) on [0, +00 [, indexed by (a, a') ExE, 
a kernel. Given two kernels F and G we define their convolution F * G as the kernel 

(i^*G'),,„.([0,t]) := V f F^^^idr) f ^ G^^^,{ds). (4.25) 

We can reduce to the case of go G {0, ... , A^}, so that Sq = 0. We set 
Z„,,,(i) :=IEa (exp(-AJ„([0,t])) l(x^^=a')) 



E„ ^exp (^^Y,vl/\n{Xk-i)^ 



l(Xiv,=«') 
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where we recall that A„(q!) := (l(j=n, o-=i) ~ l(m=i+i, o-=-i)) for = {i,cr), and 



fc=i 



Let us also set for a = («, a) and a' = {i' , a') 
Ma,a'idT) := qa,a' '^iidr) ex.p 



AA„(a) 
2r2 



qa,a' ^ exp + A A„(a)) l(^>o) dr. (4.26) 



By summing over all possible values of Nt we obtain 



(X,=a') 



fc=l 



where M^* is the convolution of M with itself £ times. Let us set for e > 

r+oo I-+00 / g 

:= e Ma,a'{dT) = qa,a' Pi V exp ( -- - (A + A A„(a)) —] dv 

= Qa,a' Ca{X, e), 

where Ca is defined in (|4.23p . and let us denote 

Ba,a' = B^ ,^, = qa,a' Pi J V exp \ -{Pi + A A„(a)) —j dv 
Because of the explicit form of the coefficients, we have 



P^ 



Pi + XAn{a) 



Bnn' = < 



{(n,+l),(n + l,-l)} 
= (n,+l) 
a = (n + 1, —1) 



a = (n,+l 



^^ /3„+i-A' 

Recall that Q := {qa,a')a,a'eE is a permutation matrix, more precisely a cyclic permutation 
of E. The matrix S is obtained by replacing two non-zero elements of Q with, respectively, 
fj^^x ^^'^ /3^"^-A • '^^^ characteristic polynomial of i? is therefore equal to 

Pn Pn+l 



p{t) = t 



2N 



Pn + \ Pn+l - A 

Recall that we assume /3„ < Pn+i- 

The case A g] — /?„,0[U]/?„+i — In this case 

PnPn+l > PnPn+1 + KPn+1 - Pn - ^) = {Pn + A)(/3„+l - A) > 0, 

and therefore B has spectral radius 

Pn Pn+l ^ ^ 



P{B) 



Pn + ^ Pn+l — A 



> 1. 
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Let US go back to the matrix B'' = {B^ ^,)a^a'<=E- In this case, all non-zero elements of B 
are modified. Indeed, if Ba,a' > 0, then B^^, = CQ,(A,e), defined as in (j4.23|) . Therefore, 
the characteristic polynomial is 

Therefore, the spectral radius of B^ is 



p{B^) = {Fn{X,e,B)) 



1 

2JV 



Since B"^ is an irreducible matrix with non-negative entries, by the Perron- Frobenius theo- 
rem [21 Th. 1.6.4], p{B'^) is an eigenvalue of B'' with multiplicity 1; moreover this eigenvalue 
is associated with a right-eigenvector (va) and a left-eigenvector (/q.) such that Tq, > and 
> for all a e E. 

Since e — > p{B^) is strictly decreasing with value p{B) > 1 at e = and limit as 
e — > -l-cxD, then there exists a unique eo such that -F„(A, eo, B) = 1. Let us set 

''a '"a 

By construction, M is a semi-Markov kernel, i.e. 

/■ + 00 

V / M„,,,(dr) = l, VaGE. 
Moreover, Za^a'it) satisfies 

oo 

Let now (X^, 'ffe-(_i)fc>o be a Markov renewal process with kernel M. Then we can write 

oo oo 

Ua,a'{[0,t]) :=Y,Mi*a'i%t]) = ^ P« (n + • • • + f, < t , 1, = a') . 
Notice that the kernel M has finite mean: 



p : = 



r+oD 

/ Tl'aMa,a'idT) e (O,0o), 



where = la'/^^^ is the unique invariant measure of X on E. Then, by the Markov 
Renewal theorem ^ Th. VII.4.3] 

hm ^f/„,«'([0,i]) = (4.27) 



Therefore we obtain 



Za,c.'{t) = Pa(ri > t) 6a,a' + e^"* Ua,a'{[0,t]) — 

r^i 



and summing over a' £ E 



E (exp (-A J„([0, t]))) = P,(ri > t) + e^«* r« ■ 
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and therefore, 

fn{\B)= lim ^ logE(exp(-AJ„([0,f]))) = eo. 
The case X G [0,/5n+i — I3n\- In this case we have p := p{B) g]0, 1]. Recall that 

oo 

= IPa(ri > t) 5a,a' + X^^^XC^'*]) ^ ^^^^^ > ^ 

where for a = {i,a) and a' = {i' , a') 

FJti >t)= (l-e~'^'^ 



and summing over a' G E 

Pi 



E (exp (-AJ„([0, t]))) > P,(ti > t) = ^1 - e"^ 

and therefore 

liminff^+oo j logE (exp (-AJ„([0, t]))) > liminft^+oo j log ^1 - e~2^^ 

~ limt^+oc 7 log (4) = 0. (4.28) 



it^+oo I lUg 1^2F^ 

Now, let us set 

Ma,a'{dT) := ■ , 

By construction, M is a semi-Markov kernel, i.e. 

hoo 



r+oo 

V / Ma,a'idT) = l, yaeE. 



a'eE 

Moreover, Za,a'{t) satisfies 

oo 

^a,a'(t) =IPa(ri > t) + ^ M^%,( [0, ) ■ 

Let now {Xk,fk+i)k>o be a Markov renewal process with kernel M. Then we can write 

oo oo 

Ua,a'i[0,t]) := X]M^>([0,t]) = Y,^aifi + --- + fe<t,X, = a'). 
1=1 e=i 

Notice that the kernel M has finite mean: 

n + OO 



r+oo 

fl:= y2 / Tl'aMa,a'(.dT) e (0,CX)), 



a,a'eE 

where u^' is the unique invariant measure of X on E. Indeed, by (j4.26p . Ma^a'idr) ~ C 
as r — > +00. Then, by the Markov Renewal theorem [21 Th. VII. 4. 3] 

hm lUc.,a'(.[0,t]) = ^. (4.29) 

Therefore we obtain 

Za,a'{t) < 6a,a' + Ua,a'{[0,t\) — 
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and summing over a' € E 

E(exp(-AJ„([0,t]))) < 1+ V Ua,a'{[0,t])^ <C{l + t), 

L — ^ T.I 



'&E 



for some constant C > 0. Hence 



fn{\B) = lim ) logE(exp(-AJ„([0,t]))) = 0. 



limsup - logE(exp(-AJ„([0,t]))) < limsup - log(C(l + t)) = 0. (4.30) 

t—f + OO t t— > + oo t 

By K28\\ and (|i30|) we obtain 

fni\B) = 

Continuity and convexity of fn{-,B). It is a standard fact that 

^. , l E{4mt])e-'Mm)) _ (E(j„([0,t])e-^^"([M)))^ 
^- IogE(exp(-AJ„([0,t]))) = (E(e-A.„(M)))^ 

is non-negative, so by passing to the limit t +00, fn{-,B) is convex and finite and 
therefore continuous. 

The Gallavotti- Cohen symmetry relation. Equation (|4.24p follows from the analogous sym- 
metry of F 

Fn{\e,B) = Fn{l3n+i-l3n-\e,B). 

□ 

Remark 4.6. The proof shows that the lack of analyticity of fn{-,B) is related to the tail 
of the distribution of or, equivalently, to the probability of having slow particles in the 
system. Indeed, the crucial estimate (j4.28p . which shows that fn > 0, follows from the 
polynomial decay of the probability that a particle takes an amount of time t to reach the 
next scatterer, namely P(rj > t) ~ t +00. This is also related to the absence of 
spectral gap of the dynamics. Physically, the origin of this phenomenon is the fact that 
the particle may get an arbitrarily small speed (taking thus an arbitrarily large amount 
of time before the next collision) which prevents the system from converging exponentially 
fast to the stationary state. This should not be regarded as an artefact of the model but 
rather as a general feature of collisional dynamics. We leave a more complete study of the 
convergence to the invariant measure to future work. 

4.4. Green-Kubo formula and fluctuation-dissipation relations. It is well-known 
that the Gallavotti-Cohen symmetry (|4.24p implies the identification of the thermal conduc- 
tivity with the variance of the time- integrated current, a relation known as the Green-Kubo 
formula. See for instance [n\ [T8] for a derivation of this fact in the context of chains of 
anharmonic oscillators, where the analyticity of the cumulant generating function is proven 
under suitable hypothesis on the interaction potentials. 

In fact, the Gallavotti-Cohen symmetry (j4.24p implies a relation between certain par- 
tial derivatives of the cumulant generating function /„, which is correct if the generating 
function is smooth for A close to 0. The identification of the partial derivatives of fn with 
the physical quantities of interest, namely the thermal conductivity and the variance of the 
current, requires an exchange of limits which is often hard to justify. 

In our models, as remarked before Proposition 14. 5 ^ the cumulant generating function is 
not smooth for A close to and Pn 7^ Pn+i and therefore the computations which are usually 
performed require some care. We are nevertherless able to obtain the desired results, by 
considering only left-derivatives when needed and by showing directly that the physical 
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quantities of interest satisfy the expected relations. Notice that we shall denote for any 
function g :] — e, e\-^ M 

g(O-) :=lim<7(t), g(0+) :=lim5(t), 
whenever any of such limit exists. 

Proposition 4.7. Let B = {Po, ^ If < l3n < /3n+i then 

dX Zn ~ 

where Ti = /3r\ see (imH . 

Proof. We consider A s] — /?n,0[. Then by Proposition 14.51 fn{X,B) > is defined by the 
relation Fn{\, fn{^,B),B) = 1. Hence by the implicit function Theorem 

^{\B) = -^(A,/„(A,e),^) / ^(A,/„(A,^),e) . 

A computation yields 
dF, 



{X,e,B)= J] C,(0,e) 



dX 

a£E\{{n,l),{n+l,-l)} 



PnPn+i / viV2- (^2 " ^i) exp (/3n + A) - - A) ^ dui dt;2, 

Jm.1 ^ \ Vl V2 2 2 y 



V a / ?;iexpf---(A + A„(a)A)^) d^; ]T C^'(A,e) 



a=(j,o-)G-B 

Since /„(0, = 0, by letting A j we find 



2 

PnPn+l /ir2 t;i 1^2 ^ {vj - vf) g-^^^-^^+i^ a!vi dv2 
= ± ^ 

E.=(,.)Gi./3jM+exp(-A^)d^ 

^ L rj. _rp 



Z^a=(j,o-)G-B Y 2 



□ 



We consider now the equilibrium case /3o = • • • = (3n = (3 > 0. 

Proposition 4.8. Let Beq :=(/?,..., /3) G /3 > 0. Then the function ]-/?,/?[ 9 A 

fn{X,Bcq) is analytic, even and 



Proof. The relation /„(A,i3cq) = fn{—^-,Bcc^ follows from the Gallavotti-Cohen symmetry 
()4.24p . The analyticity follows from the implicit function theorem. 

Since = Tn+i-, we know by Proposition 14.71 that (0,;Bcq) = 0. We consider A G 
] — /3, 0[. Then by Proposition 14.51 /„, > is defined by the relation F„(A, /n(A, ^Beq), ■Seq) = 
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1. Hence, differentiating this relation w.r.t. A and then setting A = we obtain, since 

/„(0,Seq) = ^(0,^eq)=0, 

1^(0, ^eq) = -^(0,0,5,,) / ^(0,0, Seq) . 

Now, 



(0' O'-^eq) = 0^ J^^ Vl V2 ^ (l^l - Ui)^ CXp (^-^ {vf + dvi dV2 

OF /" 1 / v'^ \ 

^(0,0,i3eq) = -2NP v-exp (^-/?yj ^^^^ = 

so that 



9A2 ' ' '^^^ /32 ivY 27r/3 iV V vr/J^' 

Let n G {0, . . . , - 2} and for < A/3 < /3, let us set 

A:=/3>0, i^{n,n + l}, /3„:=/?-^, 

and let us set B{AP) := {(3q, . . . ,f3]\f)- Notice that i3(0) = Beq- Let us set 
A:={{X,AP) : A/3G [0,/3[, A G /3 - A/?[ } 

and Qn '■ 

5n(A, A/?) := /„ (A, ^(A/3)) , V (A, A/3) G ^. 
Proposition 4.9. VFe /lave i/ie Green-Kubo relation 

''''M0,0) = -2-^%(0,0) = l.^ 



5A2 ' ' ' 9A/3 5A ' ' ' N\J ttP^' 

Proof. Notice first that gn{X,0) = /n(A,0) is analytic around A = 0, so that 

by Proposition 14. 81 Now, for fixed A/3 > 0, we have by Proposition 14.71 

^(0,A/3) = -^"'~^"+^ 
dX^'^' Zn{B) 

PnPn+l V vr V 

Therefore 



n 



dA(3 9A ' ' ' iV V 27r/35 
and the result is proven. □ 
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Remark 4.10. It is possible to prove directly that at equilibrium, i.e. = ■ ■ ■ = (^n 
/3>0, 



which shows that the formal exchange of limits in t ^ +00 and in A ^ in the formula 

^^logE(exp(-AJ„([0,t]))) 

yields indeed a correct result. Using the Gallavotti-Cohen symmetry one can prove both 
Proposition 14.91 and the further equality 

5. Confined tracers 

5.1. Generalities and physical observables. In this section, we introduce a model 
which gives rise to a qualitatively different behavior for the self-consistent temperature 
profile. The self-consistent temperature profile of the scatterers in the wandering tracers 
model was linear. We will see that in the case of confined tracers, the temperature profile 
becomes non-linear. A major difference between the two models is the dependence of the 
thermal conductivity on the set of temperatures of the scatterers. For an arbitrary temper- 
ature distribution of the scatterers, we have seen that in the case of wandering tracers, the 
conductivity was identified with a frequency of collisions of a tracer with two neighboring 
scatterers. As such and because the wandering tracer travels through the whole system, it 
was dependent on the temperature of every scatterer. In the case of confined tracers, the 
conductivity is a purely local function of the set of temperatures. 

The general structure of the process is again one of a Markov renewal process, the notations 
and proofs are strictly analogous to the case of wandering tracers. In this model, the 
disposition of the scatterers is the same but there are exactly tracer particles locked in 
between the scatterers, including the ones on the boundaries. The n-th particle moves in 
between the scatterers, in the interval In = [n— being reflected at the scatterers Wn-i 
and Wn with a random velocity p distributed according to 

4„(p)=P^/?ne-'^"'^. (5.1) 

Because the particle is reflected, the sign in the distribution is the opposite of the sign 
of the incoming velocity. Those models are described by A^ independent Markov renewal 
processes. Each scatterer exchange energy with its two adjacent tracer particles and in 
order to express the self-consistency condition, we must introduce notations to describe 
the motion of each tracer. We describe now the process describing the motion of the n-th 
particle traveling between scatterers Wn and Wn+i- The state space of the Markov chain is 
E = { — 1, +1}, with transition probability defined by qi~i = q^i^i = 0. 

Let (9n,0)Pn,o) the initial data and velocity of the n-th particle. We define Unfi = 
sign(p„_o)- We consider now the Markov chain {<Jn,k)k>o in E with initial state Xq = an,o- 
In fact, the Markov chain has a deterministic evolution cr^ fc = {—l)^o'n,o, k >0. 

For each a G E, we write a = ^{a + 1). Then the time of the first collision with a 
scatterer is 

c c ^ . n + an,o - qn,o ^ „ 

^n,0 — Jnfi{Qn,0,Pn,0) '■— > U, 

Pnfl 
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We now define the time the particle takes between two subsequent visits to the scatterers. 
Conditionally on the cr-algebra generated by {(Tn,k)k>o the sequence {Tn^k)k>i is independent 
with distribution defined by 

Pn+a,, 



(Tn,fc G or I an,k-\) = T exp 



2t 



n.k — l\-t\ 7 

1(t>0)"'^' 



(5.2) 



recall (12. 4p . The time of the A;-th collision with one of the two scatterers Wn and Wn+i is 

Sn,k '■= Sn,0 + Tn.l + • • • + r„,fc, k > 1. 

Before time 5^,0) the particle moves with uniform velocity Pn.o- Between time Sn,k-i and 
time Sn.k, the particle moves with uniform velocity and (Sn k)k>o is the sequence of 
times when qn^t £ {n,n + 1}. In particular we define the sequence of incoming velocities 
Vn,k at time Sn,k 

(Tji k 

Vn,0 ■= Pn,0, Vn,k ■= -, k > 1. 

Tn,k 

We define the stochastic process {qn,t,Pn,t)t>o with values in [n,n + 1] x M* 

{Qt,Pt) ■■= ' 



(5.3) 



(5.4) 



n + an,k-i + ^{t- Sn,k-l), ^) if Sn,k-1 < t < Sn,k, k > 1, 

Then, in analogy with Propositions 12. l l 12.21 [3TT] and 14. 1|, we have the result, 

Proposition 5.1. The process {{qn,t, Pn,t)t>o)o<n<N-i is Markov and its only invariant 
measure is given by 

N-l 



^ n ^^s^r. 



1(P„>0) fin e-^"^ + l(p„<o) e-^-+^^ 



where (3o = (3l and Pn = Pr and is the normalization constant, 



N-l 



n=0 



(5.5) 



(5.6) 



We next identify the physical quantities of interest. The energy exchanged between the 
scatterer n and its two neighboring particles during a time interval [0, t] is given by 



Eni[0,t]) :-- 



fc=0) 



fc>0:5„_i.fc<i 



recall that, by (|5.3p and (|5.4p . Vn,k and Vn,k+i are respectively the incoming and the out- 
coming velocity of the n-th particle at time Sn,k- The energy exchanged between scatterers 
n and (n + 1) during a time interval [0, t] is given by 

Jn^„+l([0,t]) := - ^ '"n.fcO-n.fc 
fc>l:5„,fe<i 

The total entropy flow S'n([0, i]) and ^([O, t]) due to the exchange of energy between the 
scatterers and a particle can be defined as in ()4.7p . The energy flow per unit time Sn in 
the stationary state, the entropy flow per unit time Sn and S, and the average current of 
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energy per unit time Jn between and 'Wn+i, can be defined as in, respectively, 
(fi:9D and (fTOjl . 

As in the case of wandering tracers (Proposition 14. 3|) , we may study the above limits 
defining the physical properties of the model. As compared to Proposition 14.31 the main 
difference resides in the expression of the energy exchanged En with the system. This is the 
origin of the difference of shapes of the temperature profiles of the two models. 

Proposition 5.2. For all n = I, . . . , N — 1, 

£n = — ^ \ ^ and So = — , £n = ^ , (5.7) 

= ^" ^"+\ 5 = V ^ " ^ > 0. (5.8) 



The proof is completely analogous to that of Proposition ()4.3p and we do not repeat it. 
The main feature of the proof is again that by using the renewal theorem, the conductivity 

_ Jn _ 1 

'^n rp rp /7 

appears as a frequency of collision of the tracer with the walls of the box to which it is 
confined. 

5.2. Self-consistency condition, temperature profile and Fourier's law. We now 

derive the consequence of the self-consistency condition iS„ = 0, n = 1, . . . , — 1 on the 
shape of the temperature profile. 
We set 

N-\ 

(x) := VIm i+ir(x)A(ri+i-ri) 
^ — ' Ljv ' jv L 



QN 

i=0 

and 



/lAr(x) :=Tl+ / gN{t) dt. 
Jo 

Notice that hN{i/N) = Ti and that Hn linearly interpolates between these values. 

Proposition 5.3. (Self-consistency condition) The only collection {Tn)n=o,...,N such 
that 

£n = 0, n = l,...,N -1 

with Tq = Tl and = Tr, is the solution of 

I {Tn-Tn+i) _^ {Tn-Tn^i) \ ^ 1 < n < A - 1. (5.9) 



V (T„)-2 + (r„+i)-2 (T„)-2 + (T„_i)- 2 
In this case, when N — > +00, /iat converges uniformly to the function 



3 



h[x):=[Tl+x{Tl-Tl)j , xG[0,l], (5.10) 
unique solution of the equation 

^ /li /i') = 0, X G]0,1[, 



^ /i(0) = Ti, h{l) = Tr 
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Proof. We first note that (|5.9p implies that at is a solution of the elliptic equation 

/ aN h'j^ ^' dx = 0, G C,°°(0,1), 

^0 



where 

Af-l 



i=0 

This can be seen by writing, 



1 1 



/ h'^ (f' dx = ^ — "T" V / f'{x)dx 



-0 ^+r,_,i 



-0 T, - + T,/, 



summing by parts this last equation and using (j5.9p . We are going to show below that the 
sequence of functions gjy is bounded in L^(0, 1). As every subsequence is also bounded, 
one can extract from every subsequence a subsubsequence, weakly converging to some 
g £ L'^{0, 1). Correspondingly, from every subsequence of h^, one can thus extract a 
subsubsequence converging uniformly (in x). Since the convergence is uniform, every limit 
function h obtained in this way satisfies 

= / aisfh'^ip'dx^ f 2h2h'ip'dx, 

i.e. 







h'^ = 0. 



with boundary conditions /i(0) = T^, /i(l) = Tr. This equation admits a unique solution 
given by (|5.10p . This implies that all the subsubsequences constructed above converge to 
(|5.10p and thus that the full sequence converges to it as well. We show now that the 
sequence gN is bounded in L^(0, 1) norm. Multiplying (|5.9p by and summing by parts, 
we obtain 

~ = -'^^i n — ^1 



_ _ 1 

n=l In J^N +^N-l ^1 +^0 

Let us suppose that Tr = Tat+i > Tq = T^. Then, because ()5.9p implies that the sign of 
the sequence (Tn+i — Tn) is constant, we have Tn+i > for all n = 0, . . . , — 1, and in 
particular Tq < T„ < T/y. One also obtains that 

_ 1 _ i_ 

rp 2 ,rp 2 

TI rp n + 1 rrp rp \ 'T 

so that the function T is increasing and concave. In particular, 

2 I rp 2 rp 2 J_ rp 2 

n=l In +J-n+l N +^N~1 



26 



RAPHAEL LEFEVERE AND LORENZO ZAMBOTTI 



and therefore Tn - Tn-i < Tl /{N -2). In particular, 

N~2 N-2 , 

n=l 

By (j5l^ we have that 

N 



1 T 2 _|_ T 2 
1 In 



< 



Tl 
Tr 



(5.12) 



»1 iV 

/ g%dx = Ny2{Tn-Tn+if <C 

•^0 n=0 



which shows that the sequence gjy is bounded in L^(0, 1). 
We define the local thermal conductivity by 



n 



lim K 



[Nx] 



for X G [0, 1]. Its numerical value is easy to obtain, 



lim K 



[Nx] 



lim Z. J I 



lim 



2 _^ f T^PlNxl+l 



h{x) 



where h is the function (|5.10p . Unlike (j4.20p . in this case the local conductivity has a spatial 
dependence. 

5.3. Cumulant generating function and the Gallavotti-Cohen symmetry relation. 

We can repeat the analysis of the cumulant generating function of subsection 14.31 in the 
context of the confined tracers. The proofs are actually simpler and the results translate 
word by word. 

6. Conclusions and prospects. 

In the tracers-scatterers models introduced in this paper, we were able to show the validity 
of Fourier's law and interpret the thermal conductivity as the collision frequency between 
tracers and scatterers. This comes naturally as a consequence of the renewal theorem 
for Markov renewal processes. We have recovered the two types of temperature profiles 
observed in deterministic systems described by local collisional dynamics. We were also 
able to study in details the cumulant generating function of the time-integrated current, 
showing in particular its lack of analyticity. We have provided a formula allowing the 
computation of cumulants of any order in and out of equilibrium. In particular we have 
shown the validity of the Green-Kubo formula in our models. 

Natural problems to study in future works are the large deviations properties and con- 
vergence to the invariant measure of those models. They may be also further extended 
and studied in different interesting ways. One possible extension is to consider cases where 
the tracers are transmitted or refiected according to some non-trivial probability ditribu- 
tion. In that case, the Markov chain associated to the Markov renewal process becomes 
non-deterministic. Although we have provided an explicit form for the invariant measure 
covering that case too, it would be interesting to study in details the dynamical properties 
of those systems. In particular, it appears that breaking the right-left symmetry in the 
transmission-refiection rules induces a modification of the stationary current which enters 
into competition with the current driven by the temperature gradient imposed in the sys- 
tem. It could be also interesting to consider models, as in [8], where the velocity of the 
tracer remains unaffected by the scatterers with some positive probability. 
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